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FUS is genetically and pathologically linked to amyotrophic lateral sclerosis (ALS) and frontotemporal lobar 
degeneration (FTLD). To clarify the RNA metabolism cascade regulated by FUS in ALS/FTLD, we 
compared the FUS-regulated transcriptome profiles in different lineages of primary cells from the central 
nervous system. The profiles of FUS-mediated gene expression and alternative splicing in motor neurons 
were similar to those of cortical neurons, but not to those in cerebellar neurons despite the similarity of 
innate transcriptome signature. The gene expression profiles in glial cells were similar to those in motor and 
cortical neurons. We identified certain neurological diseases-associated genes, including Mapt, Stxla, and 
Scn8a, among the profiles of gene expression and alternative splicing events regulated by FUS. Thus, 
FUS-regulated transcriptome profiles in each cell-type may determine cellular fate in association with 
FUS-mediated ALS/FTLD, and identified RNA targets for FUS could be therapeutic targets for ALS/FTLD. 



Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disorder characterized by selective motor 
neuron death in adulthood. The etiology of ALS remains obscure, although many pathomechanisms 
have been suggested including RNA metabolism and non-cell autonomous toxicity 1,2 . While most ALS 
cases are sporadic, approximately 10-20% of ALS patients have a family history. In this regard, several genes have 
been identified as the cause or risk-factors for ALS, such as FUS and TAR DNA-binding protein (TDP-43), that 
are RNA binding proteins known to regulate RNA splicing and transcription 3,4 . Interestingly, accumulation of 
FUS and TDP-43 is observed in the cytoplasmic inclusions in sporadic ALS and frontotemporal lobar degen- 
eration (FTLD), which is the second most common cause of presenile dementia 5,6 . Together with the evidence of 
FUS and TDP-43, which accumulate in cytoplasmic inclusions, the above studies strongly suggest that FTLD, ALS 
with dementia, and classical ALS are all parts of a clinicopathological spectrum of diseases characterized by 
disruption of RNA metabolism, which causes neurodegeneration in a subset of neurons in the central nervous 
system (CNS) 7 . 

How does loss-of- function and/or gain-of- toxicity cause neuronal cell death in diseases? Several lines of 
evidence indicate the involvement of loss-of- function in neurodegeneration observed in ALS/FTLD associated 
with FUS. Pathological studies have demonstrated FUS redistribution into the cytoplasm from the nucleus in 
sporadic and familial ALS as well as FTLD 810 . Other studies also showed the distribution of mutant FUS proteins 
associated with familial ALS in the cytoplasm from the nucleus, in contract to nuclear localization of endogenous 
and wild-type FUS U ~ 14 . Furthermore, the loss of FUS directly leads to neuronal cell death in drosophila 15 and 
zebrafish 16 . 

The cell selectivity in ALS/FTLD has remained a long-standing mystery. As mentioned above, the pathology of 
FUS-associated ALS/FTLD involves major selective neuronal vulnerability in both motor neurons and cortical 
neurons. Glial cells such as astrocytes and microglial cells are also likely to be involved in ALS/FTLD with regard 
to non-cell-autonomous toxicity 1719 . On the other hand, cerebellar neurons are typically spared in ALS/FTLD. 
The expression pattern of FUS cannot explain the cell and region specific-selectivity in ALS/FTLD, since they are 
expressed ubiquitously throughout the CNS 8 . Based on this background, it is important to define Fus-targeting 
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RNA profiles in different cell lineages from the CNS for more precise 
interpretation of selective neuronal death in FUS-associated ALS/ 
FTLD. 

We reported recently that FUS regulates a subset of exon splicing 
events and gene expression in mouse primary cortical neurons in a 
position-dependent manner 20 . We also showed scattered binding of 
FUS to and around alternatively spliced exons, including those assoc- 
iated with neurodegeneration, such as Mapt, Camk2a, and Fmrl 20 . In 
the present study, we extend our research to the global roles of FUS 
on RNA metabolism in different cell lineages from the different CNS 
regions, including primary motor neurons, cortical neurons, glial 
cells, and cerebellar neurons, and describe a new pathomechanism 
of FUS-related ALS/FTLD. 

Results 

Lentivirus-mediated silencing of FUS in different cell lineages from 
the CNS. We have described recently the establishment of Fus- 
silenced primary cortical neurons 20 . In order to compare RNA 
profiles of FUS in different cell types of the CNS, we prepared Fus- 
silenced primary motor neurons, primary cerebellar neurons, and 
primary glial cells by introducing lentivirus-expressing shRNA 
against Fus (shFUS) or scrambled control (Fig. 1A). The purity of 
each primary cell culture was confirmed by immunostaining for cell 
specific markers. We successfully established primary cultures of 
motor and cortical neurons with more than 95% purity. For 
primary glial cell cultures, we obtained astroctyes with purity of 
more than 95%, as confirmed by staining for GFAP. Both Purkinje 
cells and granule cells were the main cell populations among primary 
cerebellar neurons (Fig. IB). 

To exclude possible off-targeting effects, we used two different 
shRNAs, shFUSl and shFUS2 in experiments performed in triplic- 
ate, as described in our recent study 20 . Correlation analysis between 
shFUSl and shFUS2 showed a coefficient of determination (R 2 ) of 
0.93 in primary cortical neurons, indicating that these two shRNA 
have minimal off-targeting effects 20 . The expression levels of Fus 
mRNA were efficiently suppressed in all four types of cells by both 
shFUSl and shFUS2 by real-time quantitative PCR (Supplementary 
Fig. S1A). The FUS protein levels were decreased in all four types of 
cells by both shFUSl and shFUS2 by immunoblot (Supplementary 
Fig. SIB). Immunohistochemistry also showed markedly decreased 
FUS protein levels in all four primary cells infected with shFUSl and 
shFUS2 compared to shRNA against control scrambled oligonucleo- 
tides (shCont) (Supplementary Fig. SIC). 

Innate gene expression profiles in neurons and glial cells. The innate 
gene expression profiles of each cell-type infected with shCont were 
analyzed using the Affymetrix Mouse Exon 1.0 ST Array and 
compared by principal component analysis. The correlation coeffi- 
cients were calculated (Supplementary Table SI) and summarized in 
the correlation plot and 2D-PCA scores (Fig. 1C). We found that the 
expression profiles of cerebellar neurons are very close to those of 
cortical neurons. The expression profiles of motor neurons are also 
similar to those of cortical neurons but with a lesser extent. On the 
other hand, the expression profiles of glial cells were divergent from 
those of the other three neuronal cells. 

FUS regulates gene expressions in motor neurons, cortical neurons, 
and glial cells, but to a lesser extent in cerebellar cells. In the next 
step, we analyzed the gene expression and alternative splicing events in 
Fus- silenced and control primary cells using the Affymetrix Mouse 
Exon 1.0 ST Array. To identify the common effects of Fus-silencing in 
different cell lineages from the CNS, we compiled a list of differentially 
expressed genes and differentially spliced exons in primary motor 
neurons, cortical neurons, glial cells, and cerebellar neurons, in 
which shFUS was introduced (GEO accession numbers: GSE36153 
for cortical neurons, GSE42421 for motor neurons, glial cells, and 



cerebellar neurons). We first prepared a list of FUS-regulated genes 
which were differentially expressed by shFUS for each cell lineage by 
filtering the gene-level signal intensities with t-test p values of £0.1. 
Then, we identified differentially expressed genes shared among the 
four primary types of cells. 

The results of plot analysis of gene expression in primary cells of 
the CNS after Fus-depletion are shown in Fig. 2A. The numbers of 
differentially expressed genes with more than 1.3-fold change were 
higher in primary cortical neurons, motor neuron, and glial cells than 
in cerebellar neurons. Indeed, there were more than 2000 differenti- 
ally expressed genes in these three cell types but only 494 genes in 
primary cerebellar neurons (Supplementary Table S2). 

FUS-mediated gene expression profiles are similar between cortical 
and motor neurons and are less similar when these neuronal cells 
are compared to cerebellar neurons. We investigated the similarity 
in differential gene expression profiles regulated by FUS among 
primary motor, cortical, and cerebellar neurons. Venn diagrams 
indicated overlap in genes or exons whose expression was regulated 
in the same direction by FUS (t-test, p < 0.1). Both motor and cortical 
neurons which are destined to die in ALS/FTLD, shared appreciable 
proportions of gene expression profiles (775/2321, 33.4% of genes in 
motor neurons; 775/2470, 31.4% of genes in cortical neurons). 
Cerebellar neurons, which are considered to be spared in ALS/ 
FTLD, contained fewer numbers of FUS-regulated genes (t-test, p 
< 0.1) than motor and cortical neurons, and showed small 
overlapping rate with other neurons (58/494, 11.7% of genes with 
motor neurons; and 50/494, 10.1% of genes with cortical neurons) 
(Fig. 2B). The fold-changes in overlapped genes filtered by the t-test 
(p < 0.1) were plotted for primary motor, cortical, and cerebellar 
neurons. Comparison of gene expression profiles showed 
significantiy high correlations between motor and cortical neurons 
(R 2 > 0.80, Fig. 2C: p < 0.1, by t-test, Supplementary Fig. S2A: p < 
0.075 and 0.05, by t-test). 

FUS-mediated alternative splicing profiles are similar in cortical 
and motor neurons but largely divergent compared with cerebellar 
neurons. We next investigated the similarity in differential alternative 
splicing profiles regulated by FUS among primary motor, cortical, 
and cerebellar neurons. We filtered the exon-level signal intensities 
with t-test p values £ 0.1, and prepared a list of FUS-regulated exons 
for each cell lineage and identified shared FUS-regulated exons. In 
alternative splicing profiles, about 10% (965/10727, 9.0% of genes in 
motor neurons; and 965/8183, 11.8% of genes in cortical neurons) of 
genes showed overlap between motor and cortical neurons at the 
exon levels, whereas cerebellar neurons showed only overlap in 
about 4% (434/9907, 4.4% of genes with motor neurons; and 325/ 
9907, 3.3% of genes with cortical neurons) of genes at the exon levels 
compared with the other two neurons (Fig. 2B). The overlapped 
exons filtered by the t-test were also plotted for primary motor, cor- 
tical, and cerebellar neurons. Comparisons of exon splicing profiles 
showed a significant correlation only between motor and cortical 
neurons (R 2 > 0.7), but not between motor and cerebellar neurons, 
or cortical and cerebellar neurons (R 2 < 0.1, each) (Fig. 2C: p £ 0.1, 
by t-test, Supplementary Fig. S2B: p < 0.075 and 0.05, by t-test). 

Glial cells are similar to motor and cortical neurons in FUS-regu- 
lated gene expression profiles but not in alternative exon profiles. 

Glial cells, which are thought to be the modifier in ALS/FTLD, were 
also investigated with regard to similarity to motor and cortical 
neurons in differential gene expression and alternative splicing 
profiles regulated by FUS. Venn diagrams showed that the gene 
expression profiles of glial cells shared certain similarities with 
motor and cortical neurons (422/2074, 20.3% of genes with motor 
neurons; 453/2074, 21.8% of genes with cortical neurons, Fig. 2D). 
Comparison of gene expression profiles in the plot analysis showed 
positive correlations between glial cells and both motor and cortical 
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Figure 1 | Experimental schema and characterization of four primary cells from the central nervous system. (A) Primary motor neurons were harvested 
from the spinal cords of C57BL/6 mouse embryos at embryonic E13. Primary cortical neurons and primary glial cells were obtained from the cerebrum of 
C57BL/6 mouse at E15. Primary cerebellar neurons were obtained from the cerebellum of C57BL/6 mouse at E15. Motor and cortical neurons, but not 
cerebellar neurons, are affected in ALS/FTLD. Glial cells are disease- modifiers in ALS/FTLD. Primary cells were infected with lentivirus expressing two 
different shRNAs against FUS (shFUSl and shFUS2) and control shRNA (shCont). Total RNA was isolated and analyzed by the Affymetrix Mouse Exon 
Array. Each experiment was performed in triplicate. (B) The purity of primary motor neurons, primary cortical neurons, primary glial cells, and primary 
cerebellar neurons was confirmed by immunostaining using specific antibodies: anti-neurofilament-H antibody (SMI32R) for primary motor neurons; 
anti-R tubulin antibody (TU20) for primary cortical neurons; anti-GFAP antibody for glial cells; anti-Zicl antibody for granule cells in primary cerebellar 
neurons; and anti-calbindin antibody for Purkinje cells in primary cerebellar neurons. (C) The innate gene expression profiles of each cell type introduced 
with shCont were compared by principal component analysis (n = 3 for each cell-type). The correlation coefficients were calculated (Supplementary 
Material, Table SI) and summarized in the correlation plot (left). The 2D PCA scores and the loadings plots of the innate gene expression profiles of each 
cell-type introduced with shCont indicated that significant separation between the profiles of the three primary neurons and that of glial cells (right). 
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Figure 2 | Comparison of gene expression and exon splicing profiles among different cell lineages from the central nervous system after Fus 
knockdown. (A) Scatter plot analysis of gene expression using shCont and shFUS in different primary cells from the central nervous system. Genes whose 
expression levels changed by more than 1.3-fold after Fus-silencing with shFUSl are indicated in red. (B) The profiles of gene expression and alternative 
splicing events were compared among the three different primary neurons. Venn diagrams indicate overlap in genes ( top) and exons ( bottom) whose 
expression levels are uniquely or concordantly regulated by FUS among motor, cortical, and cerebellar neurons (p < 0.1, by t-test). (C) The fold-changes 
in overlapped genes filtered by the t-test (p < 0. 1 ) were plotted for primary motor, cortical, and cerebellar neurons. Scatter plots of fold-changes in gene 
expression levels ( top) and alternative splicing events ( bottom) after Fus knockdown. The R 2 value was calculated for genes and exons with t-test p values < 
0.1. (D) As in Fig. 2B, the profiles of gene expression and alternative splicing events in glial cells were compared with those of motor and cortical neurons. 
Venn diagrams indicate overlap in genes ( top) and exons ( bottom) whose expression levels are uniquely or concordantly regulated by FUS among motor, 
cortical, and glial cells (p < 0.1, by t-test). (E) Scatter plots of fold-changes in gene expression levels (top) and alternative splicing events (bottom) after Fus 
knockdown. The R 2 value was calculated for genes and exons with t-test p values < 0.1. 
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neurons (R 2 = 0.58 and 0.42, respectively). On the other hand, 
comparisons of exon splicing profiles did not show a significant 
correlation between motor neurons and glial cells, and cortical 
neurons and glial cells (R 2 < 0.1, each, Fig. 2E). Thus, FUS 
regulates the expression levels and alternative splicing of the largest 
number of shared genes between cortical and motor neurons. 
Similarly, though less concordantly, FUS regulates the expression 
levels of the same genes among the cortical/motor, cerebellar, and 
glial cells, but not alternative splicing of these cells. 

Characterization of genes with altered gene expression and spliced 
alternative exon regulated by FUS in different cell lineages of the 



CNS. We also analyzed Gene Ontology (GO) terms of genes shown 
in the Venn diagrams using DAVID 6.7 21-22 . The GO terms of genes 
regulated by FUS in primary motor neurons were mainly those 
involved in signaling cascades and metabolic processes that were 
similar to those in primary cortical neurons (Table 1). The GO 
terms of those in glial cells are mainly involved in the regulation of 
the immune system. The GO terms of those in cerebellar neurons 
were not available since the number of differentially expressed genes 
was too small. The GO terms of genes involved in FUS-related 
regulation of alternative splicing events in both motor and cortical 
neurons were mainly involved in various neuronal functions such as 
synapse, nerve impulse, and neuronal projection. In contrast, none of 



Table 1 | Gene Ontology terms for FUS-associated genes in each primary cell type 

GO ID Term P value 



Gene expression 
Motor neurons 



Cortical neurons 



Glial cells 



Alternative splicing 
Motor neurons 



Cortical neurons 



Glial i 



Cerebellar neurons 



GO:0008104 
GO:0006796 
GO:0006793 
GO:00 15031 
GO:0045184 
GO:0016310 
GO:0006468 

GO:00 19637 
GO:0006644 
GO:00 16055 
GO:0009100 
GO:0007264 
GO:0006650 
GO:0007242 



protein localization 
phosphate metabolic process 
phosphorus metabolic process 
protein transport 

establishment of protein localization 
phosphorylation 

protein amino acid phosphorylation 

organophosphate metabolic process 

phospholipid metabolic process 

Wnt receptor signaling pathway 

glycoprotein metabolic process 

small GTPase mediated signal transduction 

glycerophospholipid metabolic process 

intracellular signaling cascade 



0.002392 
0.002963 
0.002963 
0.0061 12 
0.006463 
0.012273 
0.014793 

3.68E-04 
4.89E-04 
5.21E-04 
5.30E-04 
5.91 E-04 
8.42E-04 
0.001227 



GO:0009615 


response to virus 


1 .68E-08 


GO:0006955 


immune response 


1.18E-07 


GO:0048525 


negative regulation of viral reproduction 


0.012033 
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defense response 


0.01839 


GO:0045087 


innate immune response 


0.026843 


GO:0008653 


lipopolysaccharide metabolic process 


0.035673 


GO:0050792 


regulation of viral reproduction 


0.047282 


GO:0006836 


neurotransmitter transport 


0.003894 


GO:0007268 


synaptic transmission 


0.018152 


GO:0032940 


secretion by cell 


0.019717 


GO:0046903 


secretion 


0.027201 


GO:00 19226 


transmission of nerve impulse 


0.028352 


GO:0007269 


neurotransmitter secretion 


0.041581 


GO:0007267 


cell-cell signaling 


0.044736 
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synapse 


6.85E-07 


GO:0042995 


cell projection 


2.54E-06 


GO:0043005 


neuron projection 


2.29E-05 


GO:0005856 


cytoskeleton 


1 .73 E-04 


GO:0005886 


plasma membrane 


1.88E-04 
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intracellular non-membrane-bounded organelle 


2.07E-04 
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non-membrane-bounded organelle 


2.07E-04 
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0.006476 
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the neuronal function terms emerged in the list of alternative splicing 
events in cerebellar neurons. Alternative splicing events in glial cells 
were mostly categorized into immunohematological functions 
(Table 1). 

Supplementary Table S3 lists the top 10 genes that were differenti- 
ally expressed in Fus-silenced motor neurons, together with the fold- 
change values in cortical neurons, glial cells, and cerebellar neurons. 
Similar to the global profile comparison shown in Fig. 2, the express- 
ion patterns of genes in the three different primary cells were similar 
except for cerebellar neurons, whose profile was less altered by Fus- 
depletion than other cell types. 

After filtering exons with genes that were differentially expressed 
by both shFUSl and shFUS2 with t-test p value of £0.1 and with 
fold-changes s 1.3-fold in each primary cell, we categorized alterna- 
tive splicing exons into category "A" to "E". In total, 44 exons were 
validated by RT-PCR and shown in Table 2, Fig. 3A-C, and 
Supplementary Fig. S3A-E. Alternative splicing exons A were spe- 
cific to primary cortical neurons and primary motor neurons; B, 



specific to primary motor neurons; C, specific to primary cortical 
neurons; D, specific to primary glial cells, E, were common among 
primary motor neurons, cortical neurons, glial cells and/or cerebellar 
neurons. 

We identified 18 alternative splicing events that were motor- and 
cortical neuron-specific, including Mapt, Dlgap4, and Snap25 
(Fig. 3A). All the validated splicing events that were motor- and 
cortical neuron-specific are shown in Supplementary Fig. S3A. We 
observed motor-neuron-specific alternative splicing events in Synjl, 
Scn8a, and Rimsl (Fig. 3B). We also identified several alternative 
spliced events that were differentially expressed in Fus-silenced cor- 
tical neurons. Fig. 3C shows representative cortical-neuron-specific 
alternative splicing events in Kcnip, Stxbpl, and Fmrl. All the vali- 
dated events are shown in Supplementary Fig. S3C. Furthermore, 
Glial-cell-specific alternative splicing events were seen in Sgce, 
Wdr35, and Fiplll (Supplementary Fig. S3C). Exons in Fxrl and 
Tsc22d2 were alternatively spliced in all primary cell types. Map4 
exonl4 inclusion was also observed in primary cortical neurons, 
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A: specific to primary cortical neurons and primary motor neurons; B, specific to primary motor neurons; C, specific to primary cortical neurons; D, specific to primary glial cells, E, common among primary 
motor neurons, cortical neurons, glial cells and/or cerebellar neurons. 



SCIENTIFIC REPORTS | 3 : 2388 | DOI: 1 0. 1 038/srep02388 



6 



r i 

■ ;v 



I - I 



tl 



I ' — I 



it. 1 1 



ii 



shCanl shIUSi Mil UV 



I [ 



iCont shFUSI SHFU52 




I 1 



u 



1 



rr 



I ■ 



I— —I 



ii 



K K 1 



11 



.li 



1 1 



IhFUSI SIlFUS! 



Kcnipl 

Ex 2 skipping 



B * * i .1 
555 



ll 



i l 



i — i 



Not Delectable in Ex2+ 



i 



Snap25 

Alternative exon (Exon 5] 

it ■ 



qPCR KM TagMan probes 



rr 



.11 



UiCont shFUSI 



gPCR with Tag Man probes 



Fmr1 

Ex 12 inclusion 




rr 



II 



il 



1 




S11FUS1 shFUS2 



Rimsl 

Alternative 3'-UTR 




I w,. 

i u. 

Not Detectable in Ex27- 




ill 


^ ' 1 1, 


I I- 
I I- 

Not Detectable 


Not Detectable in var, 203 


•Hi 


I- - h,, 


t < 




I u 

1 1- 


1 ' 

£ ahCont shFUSI iriFUS2 


a os 
8 oe 

S M 
£ 03 

§ 0 - 


i 1 1 

ihCont shFUSI INFU52 


- I 1 

fill 

1 ... 1 

o-l— ■■— — — — j- , 



Scn8a 

Alternative exon (exon 4) 



Figure 3 | Validation of representative altered splicing events in different cell lineages from the central nervous system after Fus-silencing. (A-C) 
Three representative alternative splicing events regulated by Fus are shown in each category: (A) both motor and cortical neuron-specific, (B) motor 
neuron-specific, and (C) cortical neuron-specific. All the validation of altered splicing events including category D and E, are shown in Supplementary 
Fig. S3. The top panels represent schematic splicing changes mediated by FUS. shCont and shFus resulted in splicing events shown in the top and bottom 
rows, respectively. The second panels show representative RT-PCR of the indicated exons in primary motor neurons. Similarly, the third, fourth, and 
fifth panels show representative RT-PCR of the indicated exons in primary cortical neurons, primary glial cells, and primary cerebellar neurons, 
respectively. The experiments were repeated four times using four independent sets of samples. The results of densitometric quantification of RT-PCR are 
shown in bar graphs (n = 4; mean ± SD). *p < 0.05, between shCont and shFUS (by Student's t-test). 
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motor neurons, and glial cells; but not in cerebellar neurons (Supple- 
mentary Fig. S3D). 

Interestingly, we also identified many neurological disease-assoc- 
iated genes among the profiles of gene expression and alternative 
splicing events. We list such representative genes in Table 3. 

The protein levels of representative genes with altered spliced 
events and gene expression were validated by immunoblot in all four 
primary cells in the CNS (Fig. 4). The expression of 4-repeat Tau 
(RD4) which corresponds to the exonlO ( + ) isoform of Mapt gene 
was increased by shFUS in primary cortical neurons but was unde- 
tectable in motor neurons. On the other hand, the expression of 3- 
repeat Tau (RD3) which corresponds to the exonlO (-) isoform was 
decreased in primary cortical and motor neurons. The 89kD form of 
Braf protein encoded by the exonl2 ( + ) variant of Braf gene was 
decreased in Fus-silenced primary cortical neurons. The protein 
expression level of Syntaxin-IA was upregulated in Fus-silenced 
cortical and motor neurons as observed in its mRNA levels. 

Direct binding of FUS to target mRNA was not tisseue-type speci- 
fic in CNS. To investigate the direct binding of FUS to mRNA of the 
genes and exons with altered expression, RNA immunoprecipitation 
(RIP) using differenet tissues of CNS was performed. We prepared 
the spinal cord (E13), the cerebrum (E15), and the cerebellum (E15) 
from the mouse embryos, and FUS protein was immunoprecipitated 
from each tissue (Fig. 5A). FUS-associated mRNA levels were 
evaluated by RT-PCR using specific primers for Mapt, Dlgapl, and 
Stxbpl, of which alternative splicing events were regulated by FUS. 
Primers for Gapdh and intergenic region were used as negative 
controls. The RIP results showed that the interaction levels between 
FUS and mRNA of Mapt, Dlgapl, and Stxbpl among three different 
CNS tissues were comparable, whereas these three genes showed cell- 
type specific splicing patterns in Fig 3 (Fig. 5B). There was no 
apparent binding of FUS to Gapdh mRNA or intergenic region 
(Fig. 5B). We next analyzed position dependence of FUS-binding 
to splicing targets and their effects on alternative splicing through 
comprehensive analysis of the exon array using primary glial cells and 
HITS-CLIP of mouse brain (Fig. 5C-D). We analyzed the positions of 
CLIP-tags of 121 FUS-responsive exons (29 exon skipped and 92 
exon included by shRNA) that were filtered by t-test p value S 0.1 
and fold-change of £0.67 or £1.5 for shFUSl. We combined these 
exons into a single composite pre-mRNA and prepared integrated 
RNA maps from our HITS-CLIP reads mapped to the corresponding 
genomic regions, as described in more detail previously 23 25 . The 
analysis showed that scattered FUS binding sites mainly around the 
alternatively spliced exons. Conspicuous binding of FUS was 
observed at —500 nt upstream of the 3' end of the downstream 
intron in skipped exons (arrows in Fig. 5C). This finding was 
similar to the complexity map of primary cortical neurons that was 
reported previously by our group (Fig. 5D) 20 . 



Discussion 

In the present study, we determined the gene expression profiles and 
alternative splicing events in four different primary cells of the CNS, 
with silenced Fus induced by lentivirus encoding shRNA against Fus. 
By comparing the gene expression profiles of the four primary cells, 
we found that the genes altered by Fus-silencing were fewer in cere- 
bellar neurons than in the other three cell types, although the innate 
gene expression profiles of motor, cortical and cerebellar neurons 
were similar. These findings suggest that gene expression profile 
reflects cellular response affected by Fus-silencing in each cell type, 
given that ALS/FTLD-relevant cell types showed more alterations 
than non- ALS/FTLD-relevant cell type. On the other hand, the 
profiles of alternative splicing events were only similar in motor 
and cortical neurons. Alternative splicing events in these two types 
of neurons were far different from glial cells and cerebellar neurons, 
suggesting that alternative splicing events are uniquely fine-tuned in 
a cell-specific manner (Fig. 2). Alternative splicing contributed to 
brain development in mammals during the evolution process 26,27 . 
In the adult human, the brain expresses more alternatively spliced 
transcripts than any other tissues 28 , and the majority of splicing 
events are regulated in a tissue-specific manner 29 . Our finding that 
FUS-regulated alternative splicing events were more region- and cell- 
type-specific than those of gene expression suggests that the profiles 
of alternative splicing events may reflect direct phenomenon caused 
by Fus-depletion whereas those of gene expression include both 
indirect phenomenon as well as direct suppressive effects on tran- 
scription by binding to promoter antisense transcripts 20 . Indeed, we 
observed more frequent FUS-tagged sites on genes with altered 
spliced alternative exons than genes with altered gene expression 
(Supplementary Fig. S4). In this context, the alternative splicing 
events regulated by FUS could explain the cell vulnerability in 
ALS/FTLD associated with FUS, whereas the gene expression regu- 
lated by FUS could represent the size of cellular impact in FUS- 
associated ALS/FTLD. The fact that cell death broadly occurs in 
the cerebral cortex but not in the cerebellum and that cell death 
occurs only in neurons but not in glial cells in ALS/FTLD also sup- 
ports this notion 3031 . Although we found cell-type-specific profiles of 
FUS-regulated alternative splicing events, the RIP experiments 
showed that direct interaction between FUS and target RNA was 
comparable in different CNS tissues. The complexity map showed 
similar patterns for FUS-binding positions around the alternatively 
spliced exons in cortical neurons and glial cells. These findings indi- 
cate that FUS binding to mRNA is not dependent on cell/tissue type. 
Instead, other molecules that associate with FUS in the spliceosome 
are likely to dictate cell/tissue type-specificity of FUS-mediated 
alternative splicing events. FUS protein-interaction analysis may 
provide more detail information about cell specific machinery of 
FUS on alternative RNA splicing. 



Table 3 

Gene 


List of FUS-regulated genes/exons in various neurc 
Effect by Fus-silencing Cell type 


logical disease 

Neurological Disease 


Neuronal Function 


Mapt 


Exl 0 inclusion 


motor/cortical 


FTLD 


microtuble stabilization 


Fmrl 


Exl 2 inclusion 


cortical 


Fragile X 


Translation repressor 


Fxrl 


Exl 5— 1 6 inclusion 


all 


Fragile X 


RNA-binding protein 


Stxla 


upregulation 


motor/cortical 


Williams-Beuren syndrome 


Part of the SNARE core complex 


Snap25 


Ex5 alternative exon 


motor/cortical 


Part of the SNARE core complex 


Stxbpl 


Exl 9 skipping 


cortical 


EIEE4 


binding to syntaxins 


Camk2a 


Exl 4 skipping 


cortical 


Alzheimer's disease 


LTP 


Scn8a 


Ex4 alternative exon 


motor 


EIEE13 


sodium channel 


Sgce 


Ex2 inclusion 


glial 


Myoclonus dystonia 


dystrophin-glycoprotein complex 


Fktn 


downregulation 


all 


Congenital muscular dystrophy 


glycosyltransferase 


EIEE4: epileptic encephalopathy early infantile type 4. 
EIEE13: epileptic encephalopathy early infantile type 13. 
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Figure 4 | Protein expressions of genes with altered splicing events and gene expression in different cell lineages from the central nervous system 
after Fus-silencing. Protein levels of representative genes with altered splicing events (Tau and Braf) and gene expression (Syntaxin- 1 A) were evaluated in 
different primary cells in the CNS with Fus-silencing. Expression levels of FUS, Syntaxin- 1A (Stxla), 3-repeat Tau (RD3), 4-repeat Tau (4RD), Braf, and 
actin were measured by immunoblot. 



The results showed concordant regulation of many alternative 
splicing events by FUS in both motor and cortical neurons, which 
underscores the notion that FUS-mediated ALS and FTLD could be 
considered the same disease entity based on clinicopathological and 
genetic findings 7,32 . The cell-type specific transcriptome profiles we 
established were embryonic; however many of identified FUS- 
mediated splicing events were likely to be seen in adult tissues. 
Indeed, 34% (15/44) of alternative splicing events in Table2 were 
also seen in the list of another report 33 using adult mouse brain with 
Fus-silencing (data not shown). Among them, Mapt is the most 
notable gene whose splicing was affected by Fus-depletion (Fig. 3A 
and Table 2). The inclusion exon 10 yields 4-repeat Tau (RD4), 
whereas skipping of exonlO generates 3-repeat Tau (RD3). We 
reported previously the increase of exon 10 inclusion in Mapt in 
Fus-silenced primary cortical neurons, a finding also reported by 
several other groups 20,33 35 . In this study, we also found the increase 
of exon 10 inlusion in primary motor neurons, although the protein 
level of RD4 was only dectable in cortical neurons. In this regard, it is 
intriguing that previous studies reported the presence of high RD4/ 
RD3 ratio in various neurodegenerative disorders, including 
FTLD 36,37 . However, little is known about the involvement of Tau 
pathology in motor neuron degeneration in ALS. Further studies are 
necessary to clarify the association between the FUS-Tau pathway 
and pathogenesis of ALS and FTLD. 

As described above, the entire profiles of alternative splicing 
events in motor and cortical neurons were almost identical. 
Nonetheless, we identified only motor neuron- or cortical neuron- 
specific alternative splicing events with close observation of region- 
and cell-type specific profiles of alternative splicing. Those splicing 
events may potentially represent differences in cell fate in each clin- 
ical subtype of ALS/FTLD. We identified some channel-associated 
genes, such as Synjl, Scn8a, and Rimsl as motor neuron-specific 
alternative splicing targets regulated by FUS (Fig. 3B and Table 2), 
indicating that synaptic dysfunction provoked by Fus-silencing 
seems to be one of causes of motor neuron degeneration. 

On the other hand, cortical neuron-specific alternative exons 
might affect cerebral neurons, leading to cortex pathology in 
FTLD. Stxbpl is the causative gene for epileptic encephalopathy early 
infantile type 4 (EIEE4) 38 and participates in the regulation of 
synaptic vesicle docking and fusion by associating with the 
SNARE complex, which is essential for fusion of opposing cellular 



membranes necessary for neurotransmission 39 . It is noteworthy that 
two other SNARE complex components; Syntaxin- 1 A and Snap25, 
were also identified as FUS-regulated gene and alternative exon, 
respectively (Fig. 3C and Table 2). Fmrl is essential for normal 
cognitive development and its mutation can lead to fragile X syn- 
drome characterized by mental retardation, autism, Parkinson's 
disease, and other cognitive deficits 40 . The other fragile X syn- 
drome-related gene, Fxrl, is also included among the list of alterna- 
tive exons in which Fxrl exon 16 is skipped by Fus-silencing in all cell 
types (Supplementary Fig. S3D and Table 3). 

Involvement of glial cells, such as astrocytes and microglia, likely 
modifies and exaggerates the ALS/FTLD disease process. Although it 
is not clear whether non-cell autonomous mechanism is also relevant 
to FUS-associated ALS/FTLD, we showed here that primary glial 
cells showed more altered differential genes by Fus-silencing than 
cerebellar neurons (Fig. 2A), suggesting that knock-down of FUS 
gene could have certain impact on cellular homeostasis of glial cells, 
in addition to motor and cortical neurons. 

Our global analysis also identified several genes relevant to various 
neurological diseases (Table 3). The alteration of gene expression 
and/or alternative splicing of these genes may have a large impact 
on neuronal function. It is unlikely that only one of these genes or 
exons is solely responsible for neurodegeneration in ALS/FTLD. FUS 
silencing would have a partial effect instead of total loss of function 
by altering isoforms or down/up regulation of these genes; therefore, 
it is possible that accumulation of altered genes affected by FUS- 
depletion could cause neurodegeneration after reaching a critical 
threshold level even when the individual alternative splicing event 
and gene expression are not critical by themselves. Investigation of 
FUS-targeting molecules, especially those relevant to neurological 
diseases may provide mechanistic insights into selective neuronal 
degeneration in FUS-associated ALS/FTLD. 

Methods 

Primary cells from the central nervous system. Primary motor neurons were 
harvested from the spinal cords of C57BL/6 mouse embryos at embryonic (E) day 13. 
Primary cortical neurons and primary glial cells were obtained from the cerebri of 
C57BL/6 mice at E15. Primary cerebellar neurons were obtained from cerebelli of 
C57BL/6 mice at E15. The procedure for culture of each primary cell and lentivirus 
infection are described in the online supplementary information. All animal 
experiments were performed in accordance with the National Institutes of Health 
Guide for the Care and Use of Laboratory Animals and under the approval of the 
Nagoya University Animal Experiment Committee (Nagoya, Japan). 
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Figure 5 | The binding of FUS to mRNA is not tissue and cell type specific in the CNS. (A-B) RNA immunoprecipitation (RIP) was performed to detect 
the interaction between FUS and mRNA of genes with altered splicing events using mouse cerebrum, cerebellum, and spinal cord. (A) Anti-FUS 
immunoblot of 1% input, FUS-IP, and control IP (IgG) from mouse cerebrum, cerebellum, and spinal cord at E15 embryos was shown. (B) FUS- 
associated RNA was reverse transcribed and evaluated by semi-quantative PCR with specific primers for three representative FUS-target genes with 
altered alternative splicing events, Mapt, Dlgapl, and Stxbpl. Primers for Gapdh and intergenic region (intergenic) were used as controls. (C-D) We 
compared normalized complexity map of FUS-dependent splice sites of primary glial cells and cortical neurons. shFUS-mediated alternative splicing 
events in primary glial cells (C) and primary cortical neurons (D) are compiled. Blue arrows point to conspicuous peaks at —500 nt upstream of the 3' end 
of the downstream intron. The complexity map of primary cortical neuron (D) is identical to that shown in the Supporting Information/Fig. S3 in our 
previous report 20 , which is shown for comparison with that of primary glial cells (C). 
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Lentivirus. We designed two different shRNAs against mouse Fus as well as a control 
shRNA. The targeted sequences were 5 ' -GCAACA AAGCTACGGACAA- 3 ' for 
shRNA/FUSl (shFUSl); 5'-GAGTGGAGGTTATGGTCAA-3'for shRNA/FUS2 
(shFUS2); and 5'-AATTCTCCGAACGTGTCACGT-3' for shRNA/control 
(shCont). These were cloned into a lentiviral shRNA vector (pLenti-RNAi-X2 puro 
DEST, wl6-l, a kind gift from Dr. Eric Campeau at Resverlogix Corp., Calgary, 
Alberta, Canada). Lentivirus was prepared using the protocol described by Campeau 
et al 41 . Briefly, lentiviral particles were produced in HEK293T cells by transfection 
using Lipofectamine 2000 (Invitrogen, Carlsbad, CA). The lentivirus -containing 
supernatant was collected at 48 hours after transfection, and stored at — 80°C. 
Lentivirus titer was measured using NucleoSpin RNA Virus kit (Clontech, Mountain 
View, CA). 

Microarray analysis. Total RNA was extracted from primary motor neurons, cortical 
neurons, cerebellar neurons, and glial cells by the RNeasy Mini kit {Qiagen, Hilden, 
Germany). We confirmed that the RNA integrity numbers (RIN) were all above 7.0. 
cDNA fragments were synthesized and labeled from 100 ng of total RNA using the 
GeneChip WT cDNA Synthesis Kit (Ambion, Austin, TX). Hybridization and signal 
acquisition of the GeneChip Mouse Exon 1 .0 ST exon array ( Affymetrix, Santa Clara, 
CA) were performed according to the instructions provided by the manufacturer. 
Each array experiment was performed in triplicate. The exon-level and gene-level 
signal intensities were normalized by the RMA and iterPLIER methods, respectively, 
using the Expression Console 1.1.2 (Affymetrix). We followed the gene annotation of 
the ENSEMBL version e!61, which is based on the mouse genome assembly NCBI 
build 37.1/mm 9 . All microarray data were registered in the Gene Expression Omnibus 
with accession number GSE36153 for cortical neurons, GSE42421 for motor neurons, 
glial cells, and cerebellar neurons. 

The principal component analysis of innate gene profiles in each cell type was 
conducted using the GeneSpring software {Agilent Technologies). For each cell type, 
the gene-level signal intensities of three controls treated with shCont were compared 
with those of three samples treated either with shFUSl or shFUS2, using the Student's 
t-test. The gene expression profiles for each cell-type were established with or without 
statistical filtration of the t-test p value < 0.1. We also obtained alternative splicing 
profiles by filtering the exon-level signal intensities of the probe sets on internal exons 
with t-test p value < 0.1. Comparison analysis of the profiles was completed using 
shCont and shFUSl subsets for each cell-type profile. 

RT-PCR for alternative splicing analysis. Total RNA was isolated from cells using 
RNeasy Mini Kit {Qiagen, Hilden, Germany) followed by treatment with DNasel 
(Qiagen). cDNA was synthesized from 1 ug of total RNA with the Oligo-dT primer 
(Promega, Madison, WI). The primers for each candidate exon were designed using 
Primer3 software (http://frodo.wi.mit.edu/primer3/input.htm). The sequences of 
primers are shown in Supplementary Table S4. Semi- quantitative RT-PCR was 
performed using Ex Taq (Takara Bio Inc, Otsu, Japan) at 25-30 cycles at 98 'C for 
10 sec, 60°C for 30 sec, and 72°C for 1 min. PCR products was electrophoresed on 
15% acrylamide gel and stained with ethidium bromide. The intensity of each band 
was measured by Multi Gage software (Fujifilm, Tokyo). 

Immunoblot. Cells were lysed in TNE buffer containing protease inhibitors for 
15 min on ice. The lysates were then cleared by centrifuging the cells at 13,000 g for 
15 min at 4°C. Lysates were normalized for total protein (10 jig per lane), separated 
using a 4%-20% linear gradient SDS-PAGE and electroblotted. For immunoblot, we 
used anti-FUS antibodies (A300-293A, Bethyl Laboratories, Montgomery, TX and 
4H11, Santa Cruz Biotechnology, Santa Cruz, CA), anti-Syntaxin-lA antibody 
(abeam, Cambridge, MA), anti-RD3 antibody (Millipore, Billerica, MA), anti-RD4 
antibody (Millipore), anti-Braf antibody (Thermo Scientific, South Logan, UT), and 
anti-actin antibody (Sigma, St. Louis, MO). 

RNA imunoprecipitation (RIP). Extracts were taken from mouse cerebrum, 
cerebellum, and spinal cord at El 5, normalized for total protein (1.6 mg), and applied 
for RIP using anti-FUS antibody (A300-293A, Bethyl Laboratories) and RIP-Assay 
Kit (MBL, Nagoya Japan). Immunoprecipitation using rabbit IgG was used as a 
control. Semi- quantitative RT-PCR was performed using Ex Taq (Takara Bio Inc, 
Otsu, Japan) and random primers at 22-25 cycles at 98°C for 10 sec, 60 'C for 30 sec, 
and 72 °C for 1 min. Based on the HITS-CLIP analysis information in our previous 
study, we designed primers for Mapt at exon6, Dlgapl at exonl 1, and Stxbpl at 
exon20, respectively. Primers for Gapdh and an intergenic region were used as 
controls. 

Bioinformatics analysis. The detail was described in Supplementary experimental 
procedures. 
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